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Abstract 



It has been observed earlier that, in principle, it is possible to obtain a quantum mechanical interpretation of 
higher order quantum cosmological models in the spatially homogeneous and isotropic background, if auxiliary 
variable required for the Hamiltonian formulation of the theory is chosen properly. It was suggested that for 
such a choice, it is required to get rid of all the total derivative terms from the action containing higher order 
curvature invariant terms, prior to the introduction of auxiliary variable. Here, the earlier work has been 
modified and it is shown that the action, A = P J d'^x should be supplemented by a boundary term 

in the form a — 4f3 J '^R K \fh d^x, where, '^_R, K and h are the Ricci scalar for three-space, trace of the 
extrinsic curvature and determinant of the metric on the three space respectively. The result has been tested 
in the background of homogeneous and anisotropic models and thus confirmed. 

PACS .98.80.Hw,04.50.-Fh 

1 Introduction 

Higher order theories were originally introduced to get rid of ultraviolet divergence appearing in the process of 
quantizing gravity [1]. Stell [2] has shown that Einstein- Hilbert action being modified by the introduction of 
curvature squared term in the form, J d*Xyy—g[ACijki + BR + Ci?^] is renormalizable in 4-dimensions, where, 
A,B,C are the coupling constants. Further, such a theory has also been found to be asymptotically free [3]. 
Unitarity of higher derivative quantum theory of gravity has been established by Tomboulis [4]. In an attempt 
to explain the inflationary scenario without invoking phase transition in the very early universe, Starobinsky [5] 
considered a field equation containing only geometric terms. Starobinsky and Schmidt [6] later have shown that the 
inflationary phase is an attractor in the neighbourhood of the solution of the fourth order gravity theory. Hawking 
and Luttrell [7] on the other hand observed that under suitable conformal transformation, Einstein-Hilbert action 
together with a curvature squared term turns out to be equivalent to the Einstein-Hilbert action being coupled to 
a massive scalar field. Recently, the importance of considering term in the action has further been increased, 
as it has been observed that the 4-dimcnsional Brane world [8] effective action [9] contains such term. There is 
even stronger motivation to consider higher order theory of gravity, which we discuss below. 

In analogy with the particle quantum dynamics and the quantum field theory, Hartle-Hawking [10] suggested 
that the ground state wave function of the universe should be given by 



along with the proposal that the sum should be taken over all compact four geometries. In the above expression 
hij is the metric on the three space, Ie = —i A is the Euclidean counterpart of the Einstein-Hilbert action,- 
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where, R is the Ricci scalar and 17^^ is the metric of the four-space. Tliis imphes that the universe does not have 
any boundary in the Euchdean regime. The above functional integral over all compact four geometries bounded 
by a given three geometry is interpreted as the amplitude of a given three geometry to arise from a zero three 
geometry, where ipo[hij] is finite. Thus the ground state is the amplitude of the universe to arise from nothing. 

This excellent prescription that attempted to give a probabilistic interpretation of quantum cosmology, for the 
first time, however, suffers from serious disease. The Euclidean form Ie of the Einstein-Hilbert action A is not 
positive definite and therefore the functional integral tpo[hij] corresponding to the ground state wave function of 
the universe, as proposed by Hartle-Hawking runs into serious problem, since it diverges badly. 

A remedy to the above problem was suggested by Horowitz [11]. To get a convergent integral for the grond 
state wave function ipo[hij] , he proposed an action in the form 

S=\l d^XV^[ACl^i + B{R - AKfl 

where, Cijki is the Weyl tensor, A is the cosmological constant and A, B are coupling constants. The above action 
which includes curvature squared term, is positive definite and in the weak energy limit it reduces to the Einstein- 
Hilbert action. It has there been shown that canonical quantization of the above action leads to a Schrodinger 
like equation, where an internal variable acts as a time parameter. 

Perhaps the most attractive feature of higher order theory of gravity has been established in a series of recent 
publications [12], [13], [14]. There, it has been shown that, under careful choice of auxiliary variable, required to 
cast such action in canonical form, the higher order theory of gravity, in homogeneous and isotropic minisuperspace 
model, under quantization, yields Schrodinger like Wheeler-deWitt equation with excellant features. Firstly, one 
of the true degrees of freedom of the theory disentangles from the kinetic part of the canonical variables and 
behaves as the time parameter. Next, it gives birth to an effective Hamiltonian that turns out to be hermitian 
and as a result the continuity equation identifies the nature of space and time like variables in the Robertson- 
Walker minisuperspace model under consideration. Thus, a quantum mechanical idea of probability and current 
densities in quantum cosmology is established. Further, an effective potential emerges in the process, whose 
extremization yields vacuum Einstein's equation, which is ofcourse a desirable feature in the weak energy limit of 
higher order gravity theory. In view of all the results discussed above, we understand the need for modification of 
the Einstein-Hilbert action by the introduction of higher order curvature invariant terms. 

Now, Lagrangian of a system usually involves field variables along with their first derivatives. This restriction is 
not required in classical context but somtimes becomes mandetory due to quantum mechanical requirements, such 
as positive definiteness of probability, existence of the vacuum etc. In the theory of gravitation under consideration, 
Lagrangian contains derivative of higher order, in the presence of higher order curvature invariant terms in the 
action. Thus, in such situations, in order to express the action in canonical form, that results in the Hamiltonian 
formulation of the theory, it is required to introduce auxiliary variable, in view of which the Lagrangian turns out 
to be the function of the field (along with the auxiliary one) variables and their first derivatives. For canonical 
formulation of such a theory, Ostrogradski's prescription [15] is generally followed. It suggests to fix a spacelike 
surface S in a space-time (M, 5^,^) , which is asymptotically flat. Canonical variables are the three metric hab , a 
quantity Qab realated to the extrinsic curvature tensor and their conjugate momenta, p"'^, P""^ respectively. Here, 
Qab and P"^ correspond to extra degree of freedom in higher derivative theory of gravity. The relations amongst 
<3obj-P"''j the space-time curvature R and the extrinsic curvature kab are 

0"' = 2kab 

and 

Pab = -l3Vhh"^ R 

in the absence of electric part of the Weyl tensor. Thus the trace of the extra degree of freedom Q is twice the 
trace of the extrinsic curvature and that of it's conjugate momentum P is simply related to the scalar curvature 
and the phase space is spanned by {h,Q,p, P} . Hence, Ostrogradski's prescription [15] requires as many extra 
degrees of freedoms as there are independent components of extrinsic curvature tensor. 

Boulware [16] slightly modified the above prescription and suggested a proposal that the auxiliary variable may 
be chosen as TVab, which is the momenta conjugate to the extrinsic curvature tensor kab, where, tt"^ = / P"^ di^x, 
in the absence of the electric part of the Weyl tensor and h = det{hab) ■ Both the above definitions are indeed 
confusing. It has been shown in some earlier works [12], [13], [14] that, to get a well behaved quantum description 
of higher order theory of gravity, it is mandetory to get rid of certain boundary terms from the action. Canonical 
momenta have to be defined only after doing so, and as a result, Pab in Ostrogradski's prescription and Qab in 
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Boulware's proposal will not remain proportional to the curvature scalar R. This has been explicitely worked 
out in the following section. The quantum cosmological equation corresponding to the above action, presented by 
Horowitz [11], 

9_ d d A ^ 

dhab Okab Okcd Okab 

has been constuctcd in view of the above definition of momentum. In the above Gabcd is the metric of the 
superspace that depends on hab and Aab is an artefact of the term DaDbP""^ , appearing in the Hamiltonian 
constraint equation. Now, Einstein-Hilbert action is supplemented by Gibbons-Hawking boundary term. Likewise, 
if it can be proved that action corresponding to higher order theory of gravity should be supplemented by an 
appropriate boundary term, then above quantum cosmological equation will change along with the Hamiltonian 
constraint equation and the definition of momentum. 

While constructing a minisuperspace quantum cosmological model, Horowitz [11], suggested that the auxiliary 
variable can be chosen as the 'negative of the derivative of the action with respect to the highest derivative of the 
field variable present in the action.' This prescription is again misleading, since one can then introduce auxiliary 
variable even in standard Einstein-Hilbart action, leading to completely wrong quantum dynamics. This has been 
shown in an earlier publication [12]. One can get rid of such problem if one expresses the Einstein-Hilbert action 
along with Gibbons-Hawking term. This proves the importance of boundary term in the quantum domain. The 
problem is nontrivial when Einstein-Hilbert action is replaced or modified by higher order curvature invariant 
terms in the action, since there does not exist in general, surface terms in standard form for such theories. 

Our earlier proposal [12], [13], [14] was to get rid of all removable total derivative terms from the action prior 
to the introduction of auxiliary variable. It worked nicely in the homogeneous and isotropic model. However, it 
has been observed here that, if we make a change of variable in the action by replacing the scale factor a by ^/2z 
then the total derivative terms existing in the action are different. This means total derivative terms are different 
under different choice of variables. So, our earlier suggestion of removal of all total derivative terms from the action 
prior to the introduction of auxiliary variable is not sensible enough. This implies that one has to suggest exactly 
the total derivative terms that should be eliminated from the action prior to the introduction of the auxiliary 
variable. This has been done in the present work. For an acid test of the present siiggestion we extend our work 
in some of the homogeneous-anisotropic models. It has been observed that Boulare's prescription and our earlier 
suggestion raise lot of problems in finding even the classical field equations, when applied in anisotropic models, 
that correspond to more than one field variable in the action. Question is, which of the total derivative terms 
should be removed, since one of such terms can be expressed in terms of the other. Next question is, how many 
auxiliary variables should be introduced. In the present work we have attempted to give answer to the problems 
raised above. For this purpose we have considered an action containing only curvature squared term and studied 
both the isotropic and a class of anisotropic models. 

In the following section, we have studied some problems regarding the definition of momentum that has been 
discussed above. In sections 3 and 4, wc explained the need for a boundary term to exhibit correct quantum 
description of a theory. In sections 5 and 6 we have been able to find an appropriate boundary term that would 
supplement curvature squared action, and have made a comparison with our earlier result. Sections 7 and 8 have 
been devoted to explain the problems arising in anisotropic model along with the remedy for such problems. In 
a nutshell, we have been able to prove that curvature squared action should essentially be supplemented by an 
appropriate boundary term. 



2 Problems regarding definition of momenta in i?^ Gravity 

Hamiltonian formulation of a theory requires to express the action in canonical form. It has been pointed out 
in the introduction that, for the Hamiltonian formulation of an action containing curvature squared term in the 
form, 

A = p j d'^x^f^ i^^ (1) 

the corresponding Lagrangian in a minisuperspace model involves field variables along with their first deriva- 
tives and second derivatives. So, in order to cast the action in canonical form that demands the corresponding 
Lagrangian be expressed in terms of the field variables along with their first derivatives only, one requires to 
introduce auxiliary variables. Such auxiliary variables, should not be chosen in an adhoc manner. It has been 
observed that any form of auxiliary variable produces correct classical field equations in homogeneous and isotropic 
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model. However, these have not so far been tested in anisotropic models. Further, it is not enough to get the 
correct classical field equations only, rather one should also be able to produce correct and well behaved, quantum 
description of the system under consideration. By the above statement we mean to find a variable, if it at all 
exists, with respect to which the Hamiltonian (it is the effective Hamiltonian in the present context) would turn 
out to be hermitian. With such motivation, canonical variables, as discussed in the introduction are usually chosen 
along the line of Ostrogrdski [15] as the three metric hat , and the extrinsic curvature kab on the compact three 
manifold T, along with their congugate momenta p"'^ and tt"^ respectively. Here, p°'^ contains third temporal 
derivative of the field variable and tt"*" is known to have the form tt"^ = / P"*" d^x, While, P"^ = 8f3Vh h"'' ^P, 
in the absence of the electric part of the Weyl tensor with, h = det{hab) ■ For Hamiltonian formulation of such a 
theory Boulware [16] chose tt"'' as the auxiliary variable and kab as the conjugate momentum. So, according to 
Boulware's prescription, Hamiltonian formulation of an action containing higher order curvature invariant term 
requires introduction of auxiliary variable, 

Qab ^j^ab ^^^j ^ ^ab 4^ ^3^ (2) 



which would follow from the action A, if it is defined as 

Qab ^ ab ^ ^ /gx 



To check whether these two definitions given in equations (2) and (3) are at par, we consider homogeneous and 
isotropic Robertson- Walker metric in the form 

ds2 ^ _^i2 ^ a{ty{dx^ + f{x){de^ + sin^e dcj)^)), (4) 

for which 

^P = 6(^ + ^ + A) 

k\\ = kii = k-i-i = -aa, (5) 
where, k is the curvature parameter. Now, in view of (2), 

Qll = Q22 ^ Q33 ^ 9g^2^(« + ^ + 4)« (6) 

Now the action (1) in view of (5) is, 



A = 727r'0 I {^ + ^ + ^fa'dt (7) 



It should be noted that if one varies the action with respect to the highest derivative ( o ) present in the action, as 
suggested by Horowitz [11] and followed by pollock [17], Sanyal and Modak [12] and Sanyal [13], one gets, 

- = 144.2/J(- + - + ^)a2 (8) 

which docs not resemble with the one given by equation (6) for the simple reason that the auxiliary variable here 
is the canonical momentum conjugate to the variable oe., which is different from kab ~ 

—ad. In order to obtain 
w"^ = Q"^^ through the definition given by (3) one has to choose a variable z = ^ such that 

z = ad= -kaa (9) 

The action (7), in view of this variable is expressed as, 

A = 36%/27r2/3 j {^^f zi dt (10) 

Therefore 

Q = = 1^ = 72V2 TT^/jii^ = 1447r2/?(^ + ^ + 4)a (11) 
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Thus we see that Q 's given by (6) and (11) have the same form, while they differ only by a constant factor that can 
be arranged simply by changing the pre factor to 12 instead of 8, in the definition of P"^ (and hence tt"* — Q"^ 
given in (2)). Thus, 

Qab ^ ^ab = 1213 j Vh h"'' '^R d^x (12) 

However, if one agrees with the fact that removable total derivative terms should be eliminated from the action 
prior to the definition of auxiliary variable, ie. action (1) should be supplemented by some boundary term, as 
suggested in our earlier works [12], [13], [14], then action (10) takes the following form, 

r -2 Li l2 • 
^ = 36^277^/3 / (^ + fc^ + ^)dt + 72V27r^/3fc^ . (13) 

J \/ Z Z\/ Z \/z \/ z 

In view of which, 

... 2 

Q = 7r<"* = 72^277"/?^ = 1447r2/3(- + ^)a, (14) 
s/z a 

does not match with the definition (12). Hence we find that the definition of momentum given by (2) or (12) is 
vulnarable and should not be taken seriously. 

Further, we can see that these two definitions given by (3) and either (2) or (12) do not agree in homogeneous- 

anisotropic models also. The Homogeneous and anisotropic Kantowski-Sachs, axially symmetric Bianchi-1 and 
Bianchi-111 minisupcrspacc models, can be expressed together in the following form, 

ds^ = -dt^ + a^{t)dr^ + b^{t){de^ + fl{e)d(t)'^), (15) 

where, k is the curvature index of the 2 -dimensional surface dO^ + f^{9)d(j)'^ and, 

fk = sine, 



= sinhO, 

correspond to Kantowski-Sachs metric with positive spatial curvature, k = 1, Bianchi-1, with zero spatial curvature 
k = and Bianchi-111 metric, with negative spatial curvature k = —1 respectively. Corresponding Ricci scalar 
and extrinsic curvatures are given by, 

^R = 2{- +2- + 2- + - + -). 
a b ab b-^ b'^ 

fell = -ad, k22 = kzs = -bb (16) 
Hence, In view of the definition (12), 

Tfll = Q = 4877/3- ^R, 77^2 = 77^3 = g = 4877/3o ^R (17) 

a 

So, from the prescription [16], it is clear that a pair of auxiliary variables Q and q is required to cast action (1) in 

2 

canonical form in the anisotropic model under consideration. Now under the choice of the variables x = ^ and 

.2 

y = ^ such that, 

x = aa = —kii, y = bb = —k22 = —kzz (18) 
and in view of (16), the action (1) takes the following form, 

A = 8V2wp f ^yfi+2y--^fi-yf + -]^ dt (19) 
J x y 2 X y y 

Hence, the auxiliary variables emerging from the canonical definition (3) are, 

Q = ^ = 16V2 77/?^ ^R = 1677/3- ^R 
ox \/x a 
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dA 

q= — = 32V2 -K^^ = 327r/3a (20) 
dy 

So, we observe that the auxiUary variables obtained in (17) and (20), in view of the two definitions (2) or (12) and 
(3) respectively, do not match. Further, as pointed out earlier, if one agrees that all removable total derivative 
terms should be removed from the action (1) prior to the definition of auxiliary variable, then after doing so, we 
have, 



A = sV2.t3 [yV^ii^ + 2^ + 2^(^ + 2^) - ^ - 2^ 
J X y xy X y xy^ x'^y 



7x'^ 13y^ 2x^y 2xf 



12X^ 12t/4 'ix^y 3a;y3 



Hence, 

^ 11 dA ^ r: ^ y ,X A] xil Iw^i 
= vr^i = — = 8^277/?^ - + 2^ + ^ - 

ox \/x X y xy 2y^ 

g = ;r22 = ;r33 ^ ^ 32x/27r/3v^[* +2^ + ^-1^] (22) 
oy X y xy 2x'^ 

Therefore it is noticed that different momenta and hence the auxiliary variables are revealed from the same action. 
Since, definition (3) is the canonical definition of momenta, therefore, the definition of tt"'' given by (2) or (12) 
is vulnerable and should be rejected. It should be emphasized at this stage that any attempt to quantize gravity 
starting from such a definition of momentum given by (2) or (12) is bound to produce wrong result. However, 
even after rejecting the definition (2) or (12), wc arc not at all free from ambiguity, since, definition (3) further, 
requires clarification, depending on whether action A should be supplemented by some boundary term or not, and 
if it does, it furthar depends on the specific form the boundary term. In the following toy model we show that for 
a correct and well behaved quantum dynamics, any action indeed should be supplemented by a boundary term. 



3 The impotance of boundary term in quantum theory - a toy model 

Let us take an action in the following form, 

S = ^ [^mx"^ + kx- V{x)]dt. 

It is seen that in the above action we have deliberately kept a term kx, which is essentially a removable total 
derivative term. It is not difficult to see that the classical equation of motion is not affected any way in the presence 
of such a term. However, the definition of momentum is changed, ( which is now Px = mx + k ), in view of which 
Hamiltonian changes considerabely, 

H = ^-k'-^-^ + k^ + V. 

2m. m 2 

Such an Hamiltonian leads to a Schrodinger equation that suffers from the problem of uniterity due to the presence 
of the second term which is linear in ■ This simple situation teaches us the lesson that every action should 
be supplemented by an appropriate boundary term. It has been shown earlier [12] that if we do'nt supplement 
Einstein-Hilbert action by Gibbon's-Hawking surface term then it's Hamiltonian formulation requires introduction 
of an auxiliary variable, which finally leads to wrong Wheeler-deWitt equation. Such arguments are convincing 
enough to make us believe that action containing higher-order curvature invariant terms must also be supplemented 
by a suitable boundary term. However, the requirement of such a boundary term is much apparent in the process 
of canonical quantization of R^ theory of gravity via Boulware's prescription. 



4 Canonical quantization of theory of gravity via Boulware's pre- 
scription - need for a boundary term 

Before trying to find an appropriate boundary term that should supplement an action containing curvature squared 
(i?^) term given in the form (1), let us first observe the outcome of Boulware's prescription in homogeneous and 
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isotropic model, if action (1) is not supplemented by a boundary term. As pointed out earlier, such prescription 
requires a set of phase space variables {hab,T^°''',P°''',kab}- In this prescription, hab and tt"'' are taken as the 
canonical vaiablcs whoso corresponding momenta are p"'^ and kab respectively. For such a prescription (10) is the 
action and (11) is the auxiliary variable. For the corresponding Hamiltonian formulation, let us first express the 
action (10) in the following canonical form, 

A = j [-Qz - + kQjdt + Qz, (23) 

where, M = 36\/2 tt^/?. Action (23) is canonical since the Hessian determinant S ^. is nonzero. In view of 
the above action, we observe that, 
dL 

PQ = ^ = ~^ = = ^<^b- (24) 

oQ 

Hence, Boulware's prescription is followed in the sense that, momentum, tt"** = Q has been considered to be the 
field variable, while, kab = has been taken as the conjugate momentum. It is easy to see that the same 
definition of the auxiliary variable Q is obtained from the above action along with the classical field equations. It 
is to be noted that finally we are left with a surface term Qz , that can be expressed in the following form. 



VhK^R d^x, (25) 
where, K is the trace of the extrinsic curvature. The Hamiltonian is given by, 

H = -PQ p. + - kQ (26) 

However if one is now interested in canonical quantization of the theory, it should be done with respect to the 
basic variables hab-,k°'^ instead of hab,Q- Hence, we choose z = x, which is essentially —kab, and replace pg by 
—X and Q by px • As a consequence, 

H = xp. + ^pI - kp, (27) 
The corresponding Wheeler- de- Witt equation is 

'^^z=-AM-d^^'^-xd-x ^''^ 
The interesting feature to be noted is that, here we do not require to take into account any factor ordering ambiguity, 
since x and Pz appearing in the first term on the right hand side of the Hamiltonian (27) are independent variables, 
so also z and Pz ■ Above equation can be expressed in the following form, 

= Hoi^ (29) 

where, Hq is the effective Hamiltonian. Though this equation looks like Schrodinger equation, it suffers from the 
problem that the effective Hamiltonian is not hermitian - the same problem that we have faced in the toy model 
above in section (3). Since, we have already been able to produce a hermitian effective Hamiltonian [12] for the 
same situation under consideration, therefore we claim that equation (28) or (29) do not yield the correct quantum 
description of the theory. Further, it has been pointed out earlier that unless an action is supplementad by an 
appropriate boundary term, it is always possible to introduce auxiliary variable, even in situations, where it is not 
required. In particular. It has been shown [12] that if Einstein-Hilbert action is not supplemented by Gibbons- 
Hawking term then it is possible to introduce auxiliary variable, that keeps classical field equations unchanged 
while it produces wrong Wheeler-deWitt equation. Later, the same has been shown in the context of induced 
theory of gravity [13]. This implies that boundary term plays a crucial role in quatum domain and every theory 
should be supplemented by an appropriate boundary term. Further, in the series of our earlier works [12], [13], 
it has been particularly pointed out that, if one removes all the removable total derivative terms prior to the 
introduction of auxiliary variable then the Schrodinger equation corresponding to an action containing higher 
order terms like i?^, iJ^u^i?'*'' in the homogeneous and isotropic model, exhibits some excellent features, viz., the 
effective Hamiltonian is hermitian, as a result of which quantum mechanical probabilistic interpretation of the 
theory is naturally revealed. Further, extremization of the effective potential gives classical vacuum field equation. 
In view of all these results we believe that action (1) must be supplemented by an appropriate boundary term and 
in search of that let us first remove all removable total derivative terms from action (10). 



8 



5 Removing boundary term in Boulware's prescription 

In search of an appropriate boundary term for an aetion containing curvature squared term i?,^ , given in (1), let 
us start with action (10) and get rid of aU the total derivative terms prior to the introduction of auxiliary variable. 
Thus we get action (13), 

A = M [ + k^-= + ^]dt + 2Mfc4=- 
The boundary term that supplements the above action can be expressed in the form 

a = 4l3 j VhK^Rd^x (30) 

If we again proceed as before with Q = ^ = 2M-^, given in (14), the action can be written in the following 
canonical form 

A= [ [-Qz + Mk^-^Q^ + M^]dt + Qz + a (31) 
J ^ z^/z AM s/z 

Hence the total boundary term can be expressed in the same form E, given in (25). So everything looks nice. The 
definition of the auxiliary variable and classical field equations may be realized as usual from the Euler-Lagrange 
equation, and the Hamiltonian is obtained as 

H = -,<,p,-Mjij4 + ^e^-«^^ (32, 

For canonical quantization we express the Hamiltonian in terms of the basic variables as before, for which we 
choose z = X and replace Q by and pq by —x and as a result obtain, 

iJ = ^P. + ^pJ-«J^«^-M^l. (33) 
Finally the quantum dynamical equation reads, 

y/z oz iMx ox/ zJx 

Here again we do not require to consider any factor ordering, since x and Pz appearing in the Hamiltonian (33) 

2, 

3' 



are independent variables, likewise, z and Px- With a change of variable, a=\z^ , the above equation can be 



expressed in a much convenient form, 

iTpA - _ J^w^ + kAh^ (35) 

'^da~ AMxdx^ (^)f^^'^ ^ J^' 

which can further be expressed as before, 

dip 

= Hoip 

oa 

This time, we find that the effective Hamiltonian operator Hq is hermitian and hence there exists a standard 
quantum, mechanical probability interpretation of the theory as expatiated in our earlier works [12], [13]. Further, 
we observe that we didn't have to choose an external time variable, rather an internal variable a. acts as time 
parameter. We also notice that a = ^z^ is proportianal to the proper volume , which is naturally a resonably 
good choice for the time parameter. Finally, we end up with an effective time dependent potential, 

Mk . 2 4 2 , /3Q!, 2, 

Ve= Tp^^ a;2 + fc — 3 , 

X da 2 

whose extremization leads to a power law inflationary solution of the scale factor 
« = yf(*-*o)- 

Thus, we conclude that in the early universe, if Einstein-Hilbert action is modified by the introduction of i?^ term 
then power law inflation of the Universe is a natural quantum mechanical outcome, for which one does not require 
any additional field with a potential put in by hand. Further, it appears that graceful exit might follow naturally 
as the Universe transits from quantum to classical era if Einstein-Hilbert action would have been considered in 
addition. However, this is only a proposition that has not been shown explicitely in the present work. Hence, 
we rather conclude that the action containing curvature squared term must be supplemented by a 
boundary term a, given in equation (30). Let us now compare this, with our earher work. 
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6 Comparison with our earlier result 

In our earlier works we have chosen variables different from that suggested by Boulware. Rather, we followed the 
work of Horowitz [111 in the sense that instead of Q°-^ = tt"'' = — ^4— ^ we considered, Q"'' = tt"^ = — ~ — . 

Hence wc did not require to make any change of variable and therefore, instead of (10) we started with action (7) 
and got rid of all removable total derivative terms to obtain, 

[ad'^ + - — -^^]dt + m{-d^ + 2ka), (36) 

where m = 727r^/3. It should be noted that, as a result of removing all the removable total derivative terms, the 

action now has been supplemented by a boimdary term m(|a'^ + 2ka) , which is different from a given in (30) by 
an additive factor of |md^. This is definitely awkward. A particular action must be supplemented by the same 
boundary term, independent of the choice of variable. So, there is something wrong with our earlier proposal that 
all removable total derivative terms should be removed prior to the introduction of an auxiliary variable. However, 
as already stated, if we now choose the auxiliary variable as Q = ^ = 2maa and express the action in canonical 
form, then 



(^2 _|_ l^\2 Ql 2 

-Qa + m- '- ^^\dt + Qa + mi-b? + 2ka). 

a 4ma 3 



It can easily be checked that classical field equations remain unchanged and the definition of auxiliary variable is 
recovered. However, the corresponding quantum cosmological equation changes considerably. It is now 

i7i^ = (— ^ + - — ) + — X +fc) i/j, (37) 

oa Amx ox'' x ox x 

where, a = In a and x = d, which is clearly different from that given in (35). However, this equation can again 
be expressed in the form, 

ih^ = Ho^. (38) 
oa 

The effective Hamiltonian Hq is hermitian. Further, an effective potential has been evolved whose extremization 
leads cither to vacuum Einstein's equation or inflation. This nice result that has been produced earlier [12], 
suggests that action corresponding to the curvature squared term should be supplemented by some boundary 
term m(|a^ + 2kd) which is different from a given in (30). As a result, we finally end up with a boundary 
term that contains — 967r^/3d'^ in addition to that given in equation (25). This result is not convincing. Wc 
expect that the final boundary term should be the same whether or not the action is supplemented with some 
boundary term ininitially. It is now clear that for different choice of variables we may have different removable 

2 

total derivative terms in the action. In particular we have noticed that under the choice z = ^ , the removable 
total derivative term 2md^d appearing in the action (7) gets hidden in action (10). So our earlier proposal that 
all removable total derivative terms should be eliminated from the action is clearly wrong. In addition, we observe 
that the intrinsic time parameter here, is proportional to the natural log of the scale factor, rather than the 
proper volume as obtained in section (5), and which was naturally a better choice for such parameter. These are 
mainly the two drawbacks of our earlier works. This does not mean that one can not work with the set variables 
{habiQ"'^ = 7'""'' = — 14 — }. Rather this suggests that the action (1) must be supplemented by a definite and 

appropriate boundary term. Hence, we withdraw our earlier proposal that all removable total derivative terms 
should be eliminated prior to the introduction of auxiliary variable and now suggest that the action (1) must 
be supplemented by an appropriate boundary term u, given by (30). To test our suggested prescription 
we go over to anisotropic and homogeneous cosmological models. 



7 Problem arising in anisotropic models 

Homogeneous and anisotropic Kantowski-Sachs, axially symmetric Bianchi-1 and Bianchi-111 minisuperspace mod- 
els have already been expressed altogether in equation (15), along with the expressions for Ricci scalar and extrinsic 
curvature given by equation (16). In view of the definition of Q""^ , it is clear that for the model under considera- 
tion, a pair of auxiliary variables should be introduced. If we start with Boulware's prescription and don't remove 
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total derivative terms, then the action (19) can be expressed in canonical form in view of the auxihary variables 

(20) as 

A = /[-Oi-«-^(«^ + £|)-i{i(i-ff-^)«. + i»)).i* + (39) 

The boundary term Qx + qy takes the form specified in (25) in view of the definitions of the auxiliary variables 
given in (20). Now, one can find |^ and which are supposed to reproduce the definitions of the 

auxiliary variables Q and q . Thus, 

y 1 1^ 1 y. 2 



and 



Q = 4M-^[---{-(--^)^--}], 
X 2 2 X y y 



q =16My^[^ --{-{-- ^ f --}]. (40) 
'y 4'2'x y' y" 

It shows that the definitions of the auxiliary variables Q and q given in (18) are not recovered and as a result, 
classical field equations will turn out to be wrong, although the action (39) is canonical. This implies that, Q and q 
can not be treated as independent auxiliary variables. This suggests that - only one auxiliary variable suffices 
to express the action in canonical form. Instead of (19), if one removes the removable total derivative terms 
from the action (19), and hence work with action (21), then this action can again be expressed in canonical form, 
in view of the auxiliary variables (22). However, the result is the same, ie., not the same auxiliary variables (22) 
are realized from the Euler-Lagrange equation and hence classical field equations again turn out to be wrong. 

Same result is realized if one works even with Horowitz's prescription, ie if one considers hab, hab ie. a, d and 
b,b as the basic variables instead of hab,kab- In view of the definition of given in (16), the action (1) can be 
expressed as, 

/t2"2 u 
[(— + AaP + mb + A-aVd + 26^0 + 8d66 + + 2ka + Ak%+ 

a a b b 



.a^b"^ .db^ .hb ab'^ ^ab"^ 



As it is apparent, the lagrangian corresponding to the above action (41) contains both removable and unremovable 

second derivative terms. As par our earlier prescription let us first eliminate all the removable total derivative 
terms. In the present context of anisotropic and homogeneous minisuperspace models, completely removable terms 
appearing in the action (41) are 2A;a, 4A;|6 and Ajb^b. Further, b^d can be expressed in terms of dbb and vice- 
versa. It has been observed that if one removes all these terms mentioned, and expresses rd in terms of dbb or the 
other way round and introduces two auxiliary variables Qi = and qi = — the action can be expressed 

in canonical form but the definitions of the auxiliary variables are not recovered from the Euler-lagrange equation 
resulting in wrong classical field equations. Actually it does not require to express b'^d in terms of dbb or the other 
way round as par our earlier proposal, since these are not completely removable total derivative terms. So next, if 
one removes only the three completely removable total derivative terms keeping the interchangeable terms intact 
and follow the same procedure, then it is observed that the action can not be expressed in canonical form. The 
reason has already been stated above, ie. the two auxiliary variables are not independent. 

In view of the results that have been discussed above we disregard Ostrogrdski and Boulware's prescription of 
introducing as many auxiliary variables as the number of metric co-efficients and proceed in the following section 
introducing only a single auxiliary variable. 



8 One auxiliary variable 

As, already discussed, we now introduce only one auxiliary variable to cast the action (41) in canonical form and 
to obtain the correct classical field equations. Here we follow our earlier proposal to eliminate all the removable 

total derivative terms from the action prior to the introduction of auxiliary variable. So interchangeable terms 
b-^d and dbb are kept intact and after the removal of all the three removable total derivative terms viz., 2kd, 4A;|6 
and 4|6^6 from the action (41), we are left with, 

/b^d^ ■■ ■■ b ■ 

[( + Aab"^ + Abdb -h 4-d6d + 2b^d + %dbb+ 
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4— + — + ^ + 6fc- + + 16./3[2Ma + 2-) + (42) 

It is to be noted that the boundary term 167r/3[2fc(a + 2^) + ^] is different from a given in (30) due to the 



,3 



presence of the last term viz. IGtt/?^^ . Now, let us introduce an auxiliary variable as 

dA h"^ h 

Qi = — = 167r/3[2— a + 466 + 4-a6 + 2h\ (43) 

oa a a 

In view of the above variable Qi one can express the action (42) in the following canonical form, with Mi = 167r/3, 



provided, the action (42) contains an additional term, viz., 4jrb. Notice that, this term is a total derivative 
term that was present in action (41) and so was eliminated from the action at the begining along with two other 
terms. This term contributed 167r/3^^|- to the surface term for which the surface term could not be expressed 
in the form given by a in (30). It is further noticed that, this term does not contribute to the auxiliary variable 
(43). Hence, we understand that not all total derivative terms should be eliminated from the action as proposed 
in our earlier works [12], [13], [14]. The boundary term in (42) is now Mi[2k{a + 2f)] , which can be expressed 
exactly in the form given by a in (30). So, we notice that in view of a single auxiliary variable canonization of the 
action (41) now forces us to retain required total derivative terms in the action. So, anisotropic models confirm 
our proposition that action (1) must be supplemented by a bomidary term a having the form given in (30). 

In view of this result we understand that a boundary term a should supplement the action (41), so let us write 
the correct steps to avoid confusion, if any. Action (42) should be replaced by, 



= M,J [(_ 



+ 4a6^ + 46a6 + 4-a6a + 26^a + 8dbb + A-b^b 



4— + 4— + + 6k-^ + e-^)]dt + Mi[2k{d + 2-)]]. (45) 

The auxiliary variable Qi given by (43) remains unchanged, and the action (45) can now be expressed in the 
canonical form along with the final boundary term as, 

A = M^j [-Q,{d + 2^b) + 3^Qi-'^ + ^{6b^ + k)]dt + {2kM,+Qi){d + 2^b) (46) 

The final boundary term {2kMi+Qi){d + 2^b) , appearing in (46) can be cast exactly in the form given in (25). So 
everything looks nice, the definition of Qi is recovered, so also the correct classical field equations. However, this 
choice of basic configuration space variables viz., hat, hat ie. a, d, b for which the auxiliary variable is Qi = 
(or even if we choose a, 6, b for which the auxiliary variable is gi = ^ ) put up a problem during quantization. 
The Hamiltonian in the present situation is, 

H = -[^^PaPQ. - i2M,a(l + 2fc) [^^- + - ^ - 

Now, our basic variables are (a, d, b,Pa,PcnPb)^ instead of which we worked with (a, Qi{= Pa), b,Pa,PQ,Pb)- So, 
for quantization through basic variables let us choose d = x, hence replace Qi by p^ , but since pq^ = — (d + 2|6) , 
therefore it can not be replaced simply by —x, rather it involves some more terms involving momenta Pa and pb. 
Therefore we conclude that this variable is not healthy for quantization of the theory. 

Instead of Qi = ^ , one can even choose qi = ^ , as the auxiliary variable, which in view of action (42) is 
given by, 

qi = Aba + 8ab + Sab. 

Hence, action (42) can be expressed in the following canonical form, 

. f, . ,bd g? ,6^ d6 . . , , i 

^ = / ^91(2^ + 0) - — + gi(— + 2^ + 2^) + qiindependent terms\ + Ei. 
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Not only that Ei obtained in the process is different from that given in (25), but also the auxiliary variable 
obtained from the above canonical action, in view of Euler-Lagrange equation 

qi = 46a + 8ab + 8db + 4^6^, 

b 

is clearly different from one we started with. Thus, it will again yield wrong classical field equations. This result 
is unique and astonishing, since total derivative term now directly influences classical field equations. In order to 
recover the same qi from the canonical action and as such to obtain the correct classical field equation,s one has 
to keep the total derivative term Afrb in the action, right from the begining. This is for the simple reason that 
the term 4|6^6, that did not contribute to Qi , now contributes to the auxiliary variable qi . This again proves 
our proposition that action (1) must be supplemented by a boundary term a given by (30), otherwise classical 
field equations might turn out to be wrong. The problem of quntization of the theory is not resolved with this 
variable too, for the same reason discussed above. However, in view of such wonderful result, we are now force to 
conclude that total derivative terms play a crucial role in the classical domain and not all canonical 
actions are liable to produce correct classical field equations. 

2 

It appeared in the context of isotropic model that z = ^ is a better variable, since, in view of this variable the 
removable total derivative terms that appeares in the action (13) is exactly a, given in (30). But this essentially 
due to the simplicity of isotropic model and is not true in general. In the anisotropic model with x = \,y = \, 
the removable total derivative terms in the action (21) is not simply a given in (30) but it contains some additional 
terms which is apparent from (21). However, if the action is supplemented by a given in (30) then it can be cast 
in the following canonical form in view of a single auxiliary variable Q2 = ^ , (say), and with M = Sv^tt/?, 

A = [ [-Q{x + 2x^) - Q^{^ - 3^ + 2^ + 2-) - ^]dt + Q{x + 2a;^). 
J y 2^x^ y^ xy y' AMy ' ^ y' 

The definition of Q is recovered and correct classical field equations are obtained from the Euler-Lagrange equa- 
tions. Further, the same boundary term, given in (25) is found. However, the same problem arises during 
quantization, since, pg = —x — 2x^ , therefore, quantization through basic variables can not be performed. 

At this stage, since we have been able to convince ourselves that an action containing curvature squared term 
must be supplemented by an appropriate boundary term given by a in (30) , it might appear that a pair of auxiliary 
variable might also yield correct classical field equations provided the action is supplemented by the said boundary 
term. To check this, we take action (45) and introduce a pair of auxiliary variables 

OA 6^ ■■ b ■ ■ 

Qi = — = 167r/3[2— a + 466 + 4-d6 + 2b^] 
da a a 

and 

Qi = ^ = 46a 8a6 + 8db + 4-6^. 

db 6 

Hence action (45) can now be expressed in the following canonical form, 

fr A ,■ aQ^ ^,0.^"^ ab^ , a6^ a/c^, , „ 
A = J[-Qa-qb-^ + Qi^ + -) + 6k^ + ^]dt + i:. 

In view of the above canonical action, it is readily found that same auxiliary variables are not recovered, since 
they are now, 

52 . ^ . 

Qi =4— a + 26^+4-d6 
a a 

and, 

qi = 16a6-F4^6^ -|-8d6, 
6 

which results in wrong classical field equations. Hence Ostrogradsi's prescription [15], Boulware's proposal [16] 
and Horowitz's [11] suggestions should be out and out rejected, since all of them suggested to consider the same 
number of auxiliary variables as the number of metric coefficients. Therefore, we conclude that to obtain correct 
classical field equations one has to consider only one auxiliary variable and the action must be supplemented by 
an appropriate boundary term prior to the introduction of the auxiliary variable. 
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Now, so far as the classical field equations are concern, we find that all set of variables work as well. However, 
for canonical quantization of the theory none of them is found suitable. Thus, Ostrogradski's prescription [15] 
of considering the phase space variables hab , a quantity Qab realated to the extrinsic curvature tensor and their 
conjugate momenta, p"''', P"** , for canonical quantization of higher order theory of gravity, is unacceptable in general 
and so is Boulware's prescription [16]. Hence, for the anisotropic minisuperspace models under consideration, we 
need to search for a variable such that the action expressed in terms of that variable contains a as the only 
removable total derivative term along with the possibility of having only a single auxiliary variable. 

For this purpose, let us consider a variable z = ab'^ , in view of which, the Ricci scalar (16) takes the form 

Hence, the action (1) is 

■■2 ' 



. [ .-ib r. b ^ ^Z^. ^.r.k 

A = M 4z— + Gz-^ + 4—5- - 122:-^ + 9z-j + 6kz-j + -rj-]dt + 2Mz-^, 

J z zb b^ zb^ 0'^ r r 0* 0^ 



(49) 



where, M = 167r/3. It is important to note that, now we are left with only one removable total derivative term, 
2Mi^ , can be expressed exactly in the form given by a in (30). Further there do not arise any confusion regarding 
interchangeable terms because they are now absent from the above action (49). Under such a choice of variable 
we can inaddition have only one auxiliary variable 

^ /z zb iP' , , , 

^ = ^ = 2M[--2-+3^]. (50) 

Hence, all of our requirements are thus fulfilled. The above action (49) can now be expressed in the following 
canonical form 



The final boundary term now takes the form S given in (25). In view of the above form of the canonical action 
(51), one can recover the definition of Q and find the correct classical field equations. However, the Hamiltonian 
now takes the form, 

rr ^ , b-'{bpb-2QpQ){bpb + &QPQ) ..k-'z 

H = p.PQ + z— + i2z(6^Q + 2Mfc) 'V' ^^2) 

Now under the choice i = u , we can now quantize the theory with basic variables, by replacing Q by Pu and pq 
by —u . In view of such replacements the Hamiltonian (52) can be expressed as, 

^ = -"^^ + 4M^« + 12.(6^p„ + 2Mfc) 

The Hamiltonian (53) is highly complicated due to the presence of momentum p^ in the denominator of the 
third term. Although, under quantization, z = ab"^ , which is essentially the proper volume, acts as an internal 
time parameter, as in isotropic model, it does not produce an effective Hamiltonian that is hermitian. However, 
the trouble has originated due to the presence of the two middle terms in the Ricci scalar (48). If one can find 
a variable that combines the two terms, then the complicacy of the Hamiltonian (53) will disappear, and the 
effective Hamiltonian might turn out to be hermitian. We have neither been able to find such a variable nor been 
able to prove if it exists. Nothing to worry about, even if it does not exist, since nonhcrmiticity of the effective 
Hamiltonian corresponding to curvature squared action might be a characteristic feature of any anisotropic model 
and as a result a quantum mechanical probability interpretation is not possible. It has been mentioned in the 
introduction, that unitarity of higher derivative quantum theory of gravity has been established by Tomboulis [4]. 
So, it might be conjectured along the same line that, unitarity of the effective Hamiltonian will be restored if the 
most general form of curvature squared action is taken into account. 



9 Concluding remarks 

In connection with Hamiltonian formulation of an action containing higher order curvature invariant terms, - cur- 
vature squared term in particular, in the present context, in both isotropic and anisotropic (ofcourse homogeneous) 
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minisuperspace models, a few important results have emerged, one of which might shatter our age old idea. In 
the following we list the set of results. 

Firstly, for canonical quantization of such a theory the whole superspacc is usually spanned by the phase space 
variables hab,Qab,P°''',P"'^ , according to the suggestion of Ostrogrdski [15], where 

0"*" = 2kab 



and 



Pab = -(3Vhh"'' R 



in the absence of the electric part of the Wcyl tensor, and kat is the extrinsic curvature tensor. It has been 
shown in section (2) that such a definition of momentum is not at all reliable, since it does not tally with the 
usual canonical definition of momentum. As a result, Hamiltonian constraint equation in terms of the phase space 
variables and hence the quantum dynamical equation is clearly wrong. Consideration of the electric part of the 
Weyl tensor in the definition of momentum does not change the situation, since, corresponding to the curvature 
squared action under consideration the canonical definition of momentum remains unchanged. 

Next, wc have proved the need for a boundary term for the quantum formulation of such a theory and fi- 
nally been able to show that correct Hamiltonian formulation of an action containing curvature squared term in 
particular, requires that the action should be supplemented by a boundary term a = J ^R K \fh d^x. 

Ostrogrdski's prescription [15] suggests, there should be as many auxiliary variable as the number of metric 
coefficients. It has been shown in connection with a set anisotropic models that this is not true. It is not possible 
to find correct classical field equations unless only a single auxiliary variable is chosen. Therefore, Ostrogradski's 
prescription [15] along with Boulware's proposal [16] are vulnarablc. It has been further shown that, once the 
action is supplemented by the correct boundary term a, any one of the auxiliary variables can produce correct 
classical field equation. However, if the action is not supplemented by the appropriate boundary term a, it 
might produce wrong classical field equation. This is a rather shocking result, since this goes against our age 
old conviction that total derivative terms do not contribute to classical field equations. So we end up with two 
specific conditions, viz., not more than one auxiliary variable should be taken into account and action 
containing curvature squared term must be supplemented by an appropriate boundary term a. 
Further, we understand that total derivative terms play crucial role even in classical domain and a canonical action 
do not always guarantee correct classical field equations. 

It has already been mentioned that if the two conditions are satisfied then any of the auxiliary variables is 
capable of producing correct classical field equations. However it has been shown that these variables are incapable 
of producing correct quantum dynamical equation. To obtain correct quantum dynamical equation corresponding 
to such an action, one has to find a variable such that when the action is expressed in terms of that variable, it must 
not contain any total derivative term other than a . Further, there should exist the possibility of finding only one 
auxiliary variable. This variable has been found to be the proper volume in both isotropic and anisotropic models 
that have been considered and it acts as the internal time parameter in the corresponding quantum dynamical 
equation. Hence, there exists a possibility of casting a general quantum dynamical equation corresponding to an 
action containing curvature invariant term, which will be considered in a future article. The effective Hamiltonian 
appearing in the anisotropic model suffers from the disease of nonunitarity. This might get removed when an 
action containing most general form of higher order curvature invariant terms will be considered. 

Finally, let us point out the whole scheme of canonical quantization of the curvature squared action, 

A = (3 J d^xy^ 

1. Supplement the action by a boundary term, 

a = 4:(3 j K Vh^R d^x. 

2. Replace one of the geometric variables by the proper volume z (say), and define an auxiliary variable, 

3. Express the action in canonical form in terms of the auxiliary variable Q and other field variables, and get rid 
of the excess boundary term. Thus the total boundary term is now 



= 4/3 y K Vh 



^R d^x. 
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Now, one obtains correct classical field equations and Hamiltonian formulation of the theory is thus performed. 

4. Finally, canonical quantization is performed in terms of the basic variables, which are x = z and Px along 
with other geometric field variables and their conjugate momenta. For this purpose, replace Q by Px and pq by 
X . Thus a schrodinger like equation is revealed. 

References 

W.Thirring, Phys.Rev.D77, 570 (1950). 
K.S.Stelle, Phys.Rev.D16, 953 (1977). 

E.T.Tomboulis, Phys.Lett. 97B, 77 (1950), F.S.Fradkin and A.A.Tseytlin, Nucl.Phys.B201, 469 (1982). 
E.T.Tomboulis, PRL52, 1173 (1984). 
A.A.Starobinsky, Phys.Lett.91B, 99 (1980). 

A.A.Starobinsky and H.-J. Schmidt, Class.Quant.Grav.4 695 (1987). 

S.Hawking and J.C.LuttrclL Nucl.Phys. B306, 931 (1988). 
L.Randall and R.Simdrum Pliys.Rcv.Lctt. 83, 4690 (1999). 
S.Kanno and J. Soda hcp/th-/0205188 (2002). 
J.B.Hartle and S.Hawking, Phys.Rev. D28 2960 (1993). 
G.T.Horowitz, Phys.Rev.D31, 1169 (1985). 

A.K.Sanyal and B.Modak, Phys.rev. D63, 064021 (2001), Class.Quant.Gravit, 19, 515 (2001). 
A.K.Sanyal, Phys.lett.B 542, 147 (2002). 
A.K.Sanyal, arXiv:gr-qc/0305042, Nova Science, N.Y.(2003). 
M.Ostrogradski, Mem.Acad.St.Petersbourg Series 6, 385 (1950). 

D.Boulware, A.Strominger, E.T.Tomboulis, Quantum theory of gravity, Ed. S. Christensen Adam Hilger, 
Bristol (1994). 

[17] M.D.Pollock, Nucl Phys. B306, 931 (1988). 



